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Abstract The origin of the heliospheric magnetic flux on the Sun, and hence the origin of the 
solar wind, is a topic of hot debate. While the prevailing view is that the solar wind originates 
from outside coronal streamer helmets, there also exists the suggestion that the open magnetic 
field spans a far wider region. Without the definitive measurement of the coronal magnetic 
field, it is difficult to resolve the conflict between the two scenarios without doubt. We present 
two 2-dimensional, Alfvenic-turbulence-based models of the solar corona and solar wind, one 
with and the other without a closed magnetic field region in the inner corona. The purpose of 
the latter model is to test whether it is possible to realize a picture suggested by polarimetric 
measurements of the corona using the Fe XIII 10747 A line, where open magnetic field lines 
seem to penetrate the streamer base. The boundary conditions at the coronal base are able to 
account for important observational constraints, especially those on the magnetic flux distri- 
bution. Interestingly, the two models provide similar polarized brightness (pB) distributions 
in the field of view (FOV) of SOHO/LASCO C2 and C3 coronagraphs. In particular, a dome- 
shaped feature is present in the C2 FOV even for the model without any closed magnetic field. 
Moreover, both models fit equally well the Ulysses data scaled to 1 AU. We suggest that: 1) 
The pB observations cannot be safely taken as a proxy for the magnetic field topology, as 
often implicitly assumed. 2) The Ulysses measurements, especially the one showing a nearly 
uniform distribution with heliocentric latitude of the radial magnetic field, do not rule out the 
ubiquity of open magnetic fields on the Sun. 
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1 INTRODUCTION 



Identifying the source regions of the heli ospheric magnetic flux, and h ence those of the solar wind, is a long 



standing issue in solar physics (see e.g., 



Schwenn 2006; 



Wang 



2009). The difficulties associated with this 



identification are due mainly to the difficulty of directly measuring the solar coronal magneti c field, which 
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1977 



Wang & Sheelev 



1990; 



Neugebauer et al 



1998 



2002 



Wang & Sheelev 



2006 



Abbo et al 



201C, to 



name but a few ). Although advances on polarimetric measurements w ith coronal emission lines have been 



recently made (iHabbal et al 



2001 



largely unknown (e.g., Cargill 2009) 



2003 



Lin et al 



2004 



Liu 



20091) . the coronal magnetic field remains 



Without detailed, quantitative measurements, the coronal magnetic field has been commonly con- 
structed via numerical extrapolation. While all available schemes use the photospheric magnetic field 
as boundary input, they differ substantially in how to treat the effects of electric currents on the global 
coronal magnetic field in a volume bounded by the photosphere and some outer bound ary. The elec- 



tric c u rrents may be neglected altoge ther (e.g., the potential field source surface model by 



1969; 



Schatten et al. 



Altschuler and Newkirk 



19691). they may be assu med to be purely horizontal (e.g., the current 



Zhao & Hoeksema 



Tadesse et al 



1995b . or flow exclusively along magnetic field lines 



sheet source surface model by 
(e.g., the force-free model by 
consist ently computed as a product of the plasma properties (the magnetohydrodynamic (MHD) models 



20098. or both the volumetric and sheet currents are self- 



by e.g.. lLionello et al.1120091) . Polarized brightness (pB) images of the solar corona, routinely obtained with 
space-borne and ground-based coronagraphs, often guide the extrapolation schemes such that the resul- 
tant magnetic field configuration matches available pB images. Implied here is that the density structures, 
as manifested in pB images, reflect the magnetic topology in the inner corona. And usually dome-shaped 
streamer helmets, the most prominent feature in pB images, are seen as comprising closed magnetic fields. 
It follows that the bulk of the solar wind originates from open field regions outside streamer helmets, even 



though the 



able (IKopp 



precise fraction by w hich coronal holes and the quiet Sun contribute to the solar wind is debat- 



1994 



Hu et alJ l2003. and references therein). 



However, this scenario is not universally accepted. A distinct picture has been advocated in which the 
solar wind flows along the ubiquitous open magn etic field lines that are not l imited to corona l holes or 



the quiet Sun but come from throughout the Sun (IWoo & Habbal 



1999 



2003 



Woo et al 



2004) (see also 



Woo & Druckmulleroval2008l) . Interestingly, the arguments raised to support this picture initially also came 



from density measurements. By co mbining the near-Sun pB values with radio occultation measurements 



as well as in situ solar wind data, 



Woo & Habbal 



(119991) showed that the signatures of coronal sources 
are preserved in the measured solar wind away from the Sun, contending that these density imprints are 
almost radially propagating. A further support for this scenario comes from the fact that the white light 
images of the Sun at total eclipses, when properly processed, exhibit a rich set of filamentary structures 



Woo & Habbal 



20031) . Supposing these fine 



that extend almost radially from the solar surface (see Fig. 1 in 
structures trace the magnetic field lines, this would suggest that some coronal magnetic field lines penetrate 
the dome-shaped streamer base. More importantly, these apparently open magnetic fields wer e also seen in 
the polarim etric measurements of the inner corona using Fe XIII 10747A line (see figures in 



2001 



Habbal et al 



20031) . Car efully address i ng obs ervational complications such as collisional depolarization and the 



Van Fleck effect, 



Habbal et al 



(1200 II) argued that the largely radially aligned polarization vectors may 



indeed reflect a coronal magnetic field that is predominantly radial. 

Given the importance of addressing the origin of the open magnetic flux of the heliosphere, it is sur- 
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originated by 
proposed by 



Pneuman & Kopp 



19711) (also see 



Lionello et al 



2009 



and references therein), the scenario 



Woo & Habball (120031) has not been modeled quantitatively. Here we wish to implement this 



scenario in a numerical model, thereby testing it against two fundamental constraints that the traditional sce- 
nario can readily satisfy: one is the appearance of a dome-shaped bright feature in pB images, and the other 
is the fact established by Ul ysses measurements that the radial magne tic field strength B r is nearly uniform 



logh|[l9S 



Smith et al 



refute the suggestion of a largely radially expanding solar wind ([Smith et al 



near the Sun is obviously nonuniform (e.g. 



Svalgaardet all 1978 



2001 ). Note that the latter fact was used to 



20011). as the magnetic flux 



Vasquez et al.ll20 03). Before proceeding, 



we note that from SOHO/EIT i mages, it is ob v ious th at there are a myriad of low-lying loop-like structures 



in the corona. As suggested by 



Habbal et al 



(120011) . a large-scale closed magnetic field (the "nonradial" 
component in their paper) may also help shape the large-scale corona. To simplify our treatment, we will 
simply try to answer one question: can a dome-shaped bright feature show up in a magnetic configuration 
where there is no closed field at all? In essence, this is equivalent to saying that we are interested in the 
region somewhere above the layer below which loops abound in SOHO/EIT images, and beyond which the 
contribution from closed magnetic fields is assumed to be negligible. 

In what follows, we will present two numerical models that differ in the magnetic field configuration in 
the inner corona, one with and the other without a closed field region. Both models are able to incorporate 
the essential observational constraints near the coronal base, especially the latitudinal dependence of the 
radial magnetic field. We then construct pB maps to see whether they display features similar to what is 
seen in white light observations. Moreover, model results are also compared with several crucial parameters 
observed in situ. The models are described in section|2] results from the numerical computations are given 
in section[3] and section|4]concludes this paper. 



2 DESCRIPTION OF THE NUMERICAL MODEL 

heiti Assuming azimuthal symmetry, our models solve i n the meridional plane (r, 6) the standard two 



fluid MHD equations identical to those in 



Hu et al 



1 20031) . Here r is the heliocentric distance, and 6 is the 



colatitude. The protons are heated by a flux of Alfven waves dissipated at the Kolmogorov rate Qkoi = 
p (<5tj 2 ) 3//2 /L c , where (iJu 2 ) 1 ^ 2 is the rms amplitude of velocity fluctuations of the wave field, p = nm p 
is the mass density in which n is the number density and m p the proton mass. The dissipation length L c is 

1 /2 

given by L c = g (tp) Lq (Bb tP / B) 1 where B denotes the magnetic field strength, the subscript b represents 
the value at the coronal base, Bb tP and Lq are respectively B and L c evaluated at the pole (r = 1 Rq and 
6 = 0°). In addition, g(tp) describes the dependence of L c on ip, the magnetic flux function which labels 
flow tubes. For the electrons, the classical field-aligned heat conduction is considered. While not directly 
heated, electrons receive part of the dissipated wave energy via Coulomb collisions with protons, whose 
heat conduction is neglected. With properly specified boundary conditions and given g(ip) as well as Lq, our 
numerical computation starts with an arbitrary initial state and runs until a steady state is found. This steady 
state does not rely on the initial state but is solely determined by g(ip), Lq and the boundary conditions at 
the coronal base. As a product, the computation yields a global distribution in the r — 9 plane of the proton 
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Fig. 1 Boundary conditions imposed at the coronal base, (a): Electron density n and temperature 
T. Electrons and protons are assumed to have identical temperatures at the base, (b): Wave am- 

1 /2 

plitude (<5v 2 ) and magnetic flux ip in units of the open flux rp . (c): Radial magnetic field B r 
and the total field strength B. (d): Behavior of the function g (ip) that dictates the wave dissipa- 
tion length (please see text). Model FO is given by the continuous curves while model PO by the 
dashed ones. 



function ip, the electron and proton temperatures T e and T p , as well as the wave pressure p w — p (<5u 2 ) /2. 
The magnetic field B is then obtained via B = Vip x e^/r sin 8, where represents the unit vector along 
the azimuthal direction. 



Huetal 



d2003l) . The computational domain 



The numerical scheme has been described in detail by 
extends from the coronal base (1 Rq) to 1 AU, and from the pole (colatitude 8 = 0°) to equator (8 = 90°). 
Both the pole and equator are taken as symmetrical boundaries. At the top boundary (1 AU), all dependent 
variables are linearly extrapolated for simplicity. On the other hand, at the coronal base, the proton density 
n, magnetic flux ip, temperatures (T e and T p ) and wave pressure p w are all fixed. The velocity components 
v r and vq are derived from the requirements that v is aligned with B, and mass flux is conserved along 
magnetic field lines. What distinguishes the two models lies in how the values for (n,T e ,T p ,tp,p w ) are 
imposed at the coronal base, and also in how g(ip) is prescribed. To proceed, in what follows let the two 
models be labeled FO (fully open) and PO (partially open), respectively. 

Figure[T]shows the distribution at the coronal base with colatitude 8 of (a): the density n and temperature 

1 /2 

T (the protons and electron ones are assumed to be equal at the base), (b): the wave amplitude (5w 2 ) " and 
magnetic flux given by ip/ip , ipo being the total open flux, (c): the radial magnetic field B r (thin curves) 
and magnetic field strength B (thick lines). Figured! describes the behavior of g(ip). The solid lines are 
for model FO whereas the dashed lines describe model PO. For model PO, the choice of the lat i tudina l 
profile for (n, T e , T p , ip, p w ) and that of g(ip) have been described in considerable detail in 



Hu et al 



(2003) 
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of the magnetic flux function ip. Both computations adopt the same value for the open flux ip a = 7.37 x 
10 21 Mx, whic h corresponds to an ave r age radial magneti c field of 3.37 at 1 AU, compatible with Ulysses 



measurements (iSmith & Balogh 



1995 



Smith et al 



2001). Figure indicates that in model PO the total 



flux (ip at equator) is assumed to be twice the open flux, and only the portion 9 < 60° of the solar surface 
contributes to the open flux. However, in model FO the interplanetary magnetic flux ip a is assumed to 
emerge from all over the Sun, evidenced by the fact that the magnetic flux ip at equator equals the total 
open flux ip a . Actually, the distribution of ip with 9 derives from the 9— profile of B r at the coronal base, 
presented by the thin curves in FigOJ. The specific form of ip adopted in FO reflects the observed fact that 
B r ^ oc cos 7 9 for 9 < 45° and is roughly consta nt elsewhere, establis hed by the measurements from Kitt 
Peak synoptic magnetograms at solar minima (see 



Vasquez et al. 2003 



Fig.3). 



Figure QJ indicates that the two models do not differ substantially in the profiles of n and T. While 
model PO uses a step function in T with T jumping from 1 MK to 2 MK at 9 = 60°, model FO adopts 
a somehow smoother distribution of T with a ramp connecting two values of 1 and 1.5 MK, the ramp 
spanning the range 60° < 9 < 80°. It suffices to note that, the quoted value of 1 MK is compatible with 
electron temperatures in coronal holes, and values of 1.5 and 2 MK are in line with electron temperatures in 



the quiet Sun and inside the streamer base (IHabbal et al.lll993 



Li et al.l ll998). For the n profile, model PO 



uses a distribution which is uniform in the range 60° < 9 < 90°, but decreases with decreasing 9 for the rest 
of the solar surface. The quoted values o f 1.5 x 10 8 cm" 3 at the pole and 3 x 10 8 cm~ 3 at equator both 



agree with spectroscopic measurements (IHabbal et al 



1993 



Li et al 



1998b . On the other hand, the wave 



ampli tudes given in FigQJ for both models are below the upper limit derived spectroscopically (lEsser et al. 



19991) . We note that in model PO, the waves are assumed to be absent in closed field regions. Finally, as 
there is no direct constraint on the correlation length of the turbulent Alfven waves, the g(ip) curves given 
in FigJTJi are found by a trial-and-error procedure such that the model results to be presented best match 
observations. 



3 NUMERICAL RESULTS 



The numerical results for models FO and PO are presented in the top and bottom panels of Figure [2] 
respectively. Figures and give the contours of the outflow speed v (thick black curves) superimposed 
on the distribution of the density n. Moreover, the white contours are for log n, and the background thin 
curves represent the magnetic field lines equally spaced by 0.1i/v The polarized brightness maps (in units 
of the mean disk brightness _B© of the Sun) are displayed for the heliocentric range of 2 — QRq (Figures|2j; 
and |2jl) and 4 - 30i? Q (Figures EJ: and Ef), corresponding to the Field of View (FOV) of LASCO C2 
and C3 coronagraphs on board SO HO These pB data are computed f rom t he modeled densities using 



van de Hulstl (Il950h : 



Vasquez et al 



(12003b . It is obvious from Figsf2^ 



the standard formulae given in e.g., 
and|2]3 that both models produce an outflow field characterized by a narrow wedge of slow and dense wind 
(v < 400 km/s) embedded in the fast wind. In model PO, a region associated with closed magnetic field can 
be seen whose tip is located at around 3.5 Rq. However, model FO does not possess such a region, the solar 
wind flows from throughout the Sun. Despite this, both models produce similar pB maps in the two different 
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Fig. 2 Global solar wind models where the magnetic field is fully open (top panels, model FO) 
and partially open (bottom panels, model PO). (a) and (b): Outflow speed v (thick black con- 
tours) superimposed on the density map (logarithm of density log n). The white curves represent 
the log n contours. The Background thin black curves are magnetic field lines equally spaced by 
O.I^q, tp being the open flux. The remaining panels give the polarized brightness (pB) maps 
computed from the modeled densities n. pB is measured in units of Bq, the solar mean disk 
brightness. Note that pB(r / Rq) 2 is plotted instead of pB. Panels (c) and (d) are for the helio- 
centric range of 1— 6Rq, while (e) and (f) are for 4—30 Rq. Note that the two regions correspond 
to the Field of View for LASCO C2 and C3, respectively. 



PO this feature is associated with the trapping of plasmas in closed field regions, in model FO it is due to 
the fact that the density in the near-equator region decreases more slowly with radial distance than n does in 
the polar region. Moreover, the high-latitude magnetic filed lines can be seen to bend towards the equator to 
relax the latitudinal gradient in the magnetic pressure. From this we conclude that, even though the magnetic 
field lines in model FO are not purely radial, the surprising linear polarization measur ements with Fe XIII 
10747A line placed in the context of pB measurements (e.g., the west limb in Fig.l of 



Habbaletal 



2001) 



can be understood: the large-scale bright features may not reflect closed field regions but are associated with 
open field lines that penetrate them. 

Figure|3]compares the latitudinal distribution at 1 AU of (a) the outflow speed v, (b) flux density nv and 
(c) the radial magnetic field B r , derived from the two models. Model PO (FO) is described by the dashed 
(solid) curves. Also plotted are the daily averages of Ulysses data during the first half of the fast latitudinal 
scan from Sep 12 1994 to Mar 4 1995. The two models, although having distinct magnetic configurations 
in the inner corona, provide equally good fits to the Ulysses data. One may say that the proton flux density 
nv in model FO seems to be slightly larger than th e measured values ne ar the equator, it is nonetheless 



within the ranges typically measured in situ (Fig. 4 in lMcComas et alJl2000h . Concerning the distribution of 



B r , model FO actually agrees with the Ulysses measurements better than model PO does. Furthermore, it 
is interesting to see that B r hardly displays any latitudinal dependence in the fast wind region. But this is 



Magnetic Topology of Solar Corona 



7 



1000 






BOO 


(a) 




600 




"V t\\| 
V 1 \ \ - 


400 






200 














6 
5 

C 4 



I 2 




20 40 

en 



20 40 




20 40 60 

en 



Fig. 3 Comparison with Ulysses data of derived latitudinal distributions of various parameters, 
(a): Outflow speed v. (b): Flux density nv. (c): Radial component of the magnetic field B r . Model 
FO is given by the continuous lines, whereas model PO by the dashed lines. Moreover, in situ 
measurements by Ulysses are given by the light continuous lines. 



can be understood from the force balance condition in the direction perpendicular to the ambient magnetic 
field. When the magnetic field becomes largely radial, this condition reduces to that (/? + 1) B 2 /8ir is 
latitude independent, where (3 = 8it (p e + p p + p w ) /B 2 , p e and p p are respectively the electron and proton 
pressures, and p w the wave pressure. It turns out that this happens at some distance where /3 < 1 holds in 
both models. Consequently, the magnetic field should exhibit little latitudinal variation from there on. This 
relaxation process necessarily leads to a significant non-radial component of the solar coronal magnetic 
field, in view of the base distribution of B r which varies significantly from the pole to equator. However, 
a uniform latitudinal distribution of B r at 1 AU does not necessarily invalidate a scenario where open 
magnetic fluxes originate from throughout the Sun. In this regard, we agree with 
concluded that the coronal magnetic field cannot be exactly radial. 



Smith et al 



(120011) who 



4 CONCLUDING REMARKS 

Identifying the origin of the Sun's open magnetic flux is of crucial importance in establishing the connection 
of the in situ solar winds to their sources. In the absence of a definitive measurement of the solar coronal 
magn etic field, this identification problem is subject to considerable debate (e.g 



Schwenn 



2006 



Wang 



20091) . This is true even at solar minima when the configuration of the solar corona is relatively simple, the 
most prominent feature being the bright streamer helmets in white light images. While the prevailing view 



is that the majority of the solar wind originates from outside streamer helmets (IPneuman & Kopp 



1971). 



there also exists the suggestion th at the open magnetic field is ubiquitous on the Sun and not confined 



to coronal holes or the quiet Sun (IWoo & Habbal 



numerical model has been a common practice (e.g. 



1999 



2003 ). Imple menting the former scenario in a 



Lionello et al 



2009), however, so far the latter has not 



been modeled quantitatively and hence tested quantitatively. Here we offer such an implementation. 

We have constructed two 2-dimensional, Alfvenic-turbulence-based models of the solar corona and 
solar wind, one with and the other without a closed magnetic field region in the inner corona. The purpose 
of the latter model is to minimize the contribution of the closed magnetic field, trying to mimic a corona 
permeated with open magnetic fields which may infiltrate the dome-shaped streamer base. In specifying 
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constraints, especially those on the magnetic flux distributions. Interestingly, the two models provide similar 
polarized brightness (pB) maps in the field of view (FOV) of the SOHO/LASCO C2 and C3 coronagraphs. 
In particular, a dome-shaped feature is present in the C2 FOV even for the model without any closed 
magnetic field. Moreover, both models fit equally well the Ulysses data scaled to 1 AU. Hence we suggest 
that: 1) The pB observations cannot be safely taken as a proxy for the magnetic field topology, as often 
implicitly assumed. 2) The Ulysses measurements, especially the one indicating that the radial magnetic 
field strength is nearly uniformly distributed with heliocentric latitude, do not rule out the ubiquity of open 
magnetic fields on the Sun. 

We do not intend to resolve the conflict of the two distinct scenarios currently available for the origin 
of the heliospheric magnetic flux. Rather, the presented numerical results suggest the likelihood that the 
magnetic field structure of bright features (e.g., helmet streamers) in the corona may be more diverse than 
traditionally viewed: the magnetic flux therein can be either closed or open. To differentiate the scenar- 
ios, it is likely that more stringent constraints come from the SOHO/UVCS measurements. For instance, 
measurements based on the Doppler dimming technique have yielded that along the direction transverse to 
the streamer helmet, there exists a transition in the inferred plasma speed from unmeasurable to sig nificant 



values, and this transition seems to trace the streamer legs dStrachan et al . 2002; 



Frazin et al. 2003]), identi- 



Kohl etal 



1997). Therefore, it 



fied by the enhancement of the intensity ratio of O VI A 1032 A to H I Lya 
remains to be seen whether a model permeated with open magnetic fluxes can account for this feature. To 
do this, an obvious need is to incorporate O s+ ions in a three-fluid model and test both sc enarios. At the 



moment, in such m odels only the traditional partially open scenario has been adopted (e.g., 



Li et al 



2006 



Ofman et alj|201 ll) . Instead of implementing a further construction, let us end here 



also ask whether these observational features (IKohl et al 
universal for all streamers. 



1997 



Strachan et al.ll2 002: 



by noting that one 



may 



Frazin et al. 2003) are 
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